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Designing Small Silicon Quantum Dots with Low Reorganization Energy 
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A first principles, excited state analysis is carried out to identify ways of producing silicon quantum 
dots with low excitonic reorganization energy. These focus on the general strategy of either reducing 
or constraining exciton-phonon coupling, and four approaches are explored. The results can be 
implemented in quantum dot solids to mitigate polaronic effects and increase the lifetime of coherent 
excitonic superpositions. It is demonstrated that such designs can also be used to alter the shape of 
the spectral density for reorganization so as to reduce the rates of both decoherence and dissipation. 
The results suggest that it may be possible to design quantum dot solids that support partially 
coherent exciton transport. 


I. INTRODUCTION 

Recent progress in understanding the electronic wave 
functions of colloidal semiconductor quantum dots (QDs) 
has led to an ability to synthesize isolated, quantum 
confined building blocks with a variety of tailored op¬ 
tical properties. Quantum dots composed of silicon 
(SiQDs) seem particularly promising going forward for 
many biomedical, display, computing and solar energy 
applications. This is because they are environmen¬ 
tally benign, resource plentiful and benefit from decades 
of industrial know-how. SiQDs have been shown to 
carry out multiple-exciton generation (MEG),^’^^’^^ ef¬ 
ficiently transport excitons,^^ have higher coupling for 
same surface-to-surface separation,are free of defects^^ 
and resist oxidation betterSiQDs encapsulated within 
an inorganic amorphous matrix^^~^^ and within an or¬ 
ganic polymer blend^^“^^ for photovoltaic application has 
been studied. 

For many applications, such as light emission,^ opti¬ 
cal computing,^, and biomedical imaging, sensing and 
treat ments,^'^ carrier transport is not a primary con¬ 
cern. Within the solar energy arena, though, quantum 
dots solids need to exhibit efficient carrier dynamics. 
Transport of energy and charge is now a central issue in 
realizing the mesoscopic potential of quantum dot solids. 
This is a critical bottleneck because charge and exciton 
transport tend to proceed via low mobility, incoherent 
hopping associated with polaronic trapping and weak 
electronic coupling. 

For quantum dot solids in which excitons readily dis¬ 
sociate, a number of promising strategies seek to improve 
charge transport by focusing on properties such as trans¬ 
lational symmetry, electronic overlap, matrix encapsula¬ 
tion, and crystalline orientation. Dynamics within mate¬ 
rials composed of small SiQDs, though, tend to be domi¬ 
nated by the motion of excitons though. While quantum 
confinement in SiQDs offers benefits, such as a pseudo- 
direct gap and high emission cross section, it also results 
in a low dielectric constant. This implies low screening 
and a high excitonic binding energy so that assemblies of 
SiQDs support excitons with a Frenkel character. Their 
dissociation is sufficiently problematic to motivate a con¬ 
sideration of ways in which the excitons themselves can 


be used to efficiently conduct energy. This is reminis¬ 
cent of exciton transport that is efficiently carried out in 
photosynthetic complexes, offering either biomimetic or 
bio-inspired solution strategies. 

In an ideal assembly of such dots, photon absorption 
results in a coherent superposition of Frenkel excitons 
that will subsequently exhibit ballistic travel through 
the quantum dot solid. In reality, though, such coher¬ 
ence and the associated wavelike transport is rapidly lost. 
This is because the photoexcitation of a quantum dot 
shifts its equilibrium geometry and the dot shape will re¬ 
organize into a new configuration, as shown in Figure 1. 
The energy change associated with such a relaxation is 
a measure of exciton-phonon entanglement known as the 
reorganization energy, A. Such exciton-photon entangle¬ 
ment destroys the quantum phase coherence among the 
excitonic superposition, and transport reduces to a dif¬ 
fusive random walk of excitonic states with a statistical 
distribution of occupation probabilities. 


Ground State Excited State 



FIG. 1. Reorganization following QD excitation for SissHae. 
Ground state structure (left) relaxes into a new structure 
(right) when promoted into one of its degenerate first excited 
singlet states. As a result, the dot contracts along the vertical 
axis. The displacement is magnified 5 times in order to make 
this distortion clear. 

In the absence of any structural disorder, this reduc¬ 
tion to diffusive motion results in a slower rate of ex¬ 
citonic transport. The coherent wave packet, though, 
can quickly become trapped in naturally occurring disor¬ 
dered regions^^, so that the reduction to diffusive trans¬ 
port may improve the overall transport rate in realistic 
QD solids. 
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Intriguingly, it may be possible to combine the best 
of both transport regimes so that a superposition of 
Frenkel excitons moves in a partially coherent man¬ 
ner. Inspiration for such a goal comes from photo¬ 
synthetic complexes^^’^^ wherein proteins and pigments 
are exquisitely structured so the rate of phonon en¬ 
tanglement is slow and the hopping transport is high- 
prerequisites for exciton transport with a degree of wave- 
like character^^^^^ . At issue is whether or not such a 
balance can be achieved in QD solids. 

First consider the Hamiltonian of an isolated quantum 
dot, restricted to be in either its ground state, referenced 
to zero energy, or first excited singlet state with energy, 
Sq. If is the dimensional (length x y^mass) amplitude 
operator of vibrational mode ^ with associated momen¬ 
tum operator P^, their non-dimensional counterparts are 


given by = 
nian is then 


X and h 




P^. The Hamilto- 


^(dot) ^ J2’^{p‘l+ccUl)+c^c{eo+J2 

( 1 ) 

where c is the fermionic exciton annihilation operator and 
the equilibrium amplitude operator of each vibrational 
mode is qo^. This can be re-grouped and simplified, using 
the bosonic phonon annihilation operator, + 

ip^)/\/2, into components associated with excitonic site 
energy, Hx, phonon energy, Pph, reorganization energy, 
Preorg, ^ud excitou-photou coupling energy, Px-ph- 

^(dot) =Hx+ 4h + ^reorg + Hx-ph (2) 

where 

Hx = eoc^c 

= + (3) 

Preorg — ^ ^ > huJ^qQ^C^C (4) 

^x-ph = --^c^c^hw^qo^(bl + b^) 

The reorganization energy (Figure 2) in the position ba¬ 
sis is therefore given by A := 

Huang-Rhys factors are obtained by dividing its spec¬ 
trally resolved components, A^, by huj^. This additive 
decomposition lends itself to a continuum description of 
the exciton-phonon coupling as a smooth function of vi¬ 
brational frequency. 


3{uj) = 27r huj^X^6{uj — cj^), (5) 


known as the spectral density function. A broadened ver¬ 
sion of this function is approximated as a superposition 



FIG. 2. (Left) Vertical excitation to an excited state mani¬ 
fold, due to photon absorption, is followed by reorganization 
of the nuclei. The overall displacement can be spectrally re¬ 
solved into that of individual phonon mode amplitudes, q^. 
The excited state manifold is assumed to be a displaced ver¬ 
sion of its ground state counterpart.(Right) Excited state sur¬ 
faces depend on spin state as depicted for the pair of orbitals 
comprising an exciton. 


of shifted Drude-Lorentz distributions:^^’^^ 


N 




2 A^y^^ 




(6) 


where and Qi characterize the width and shift of the 
^th L)pxide-Lorentz peak respectively. The rate at which 
structural reorganization occurs following excitation can, 
in large measure, be anticipated by identifying two fea¬ 
tures of this distribution. This will be considered after 
constructing a Hamiltonian that allows for exciton trans¬ 
port through an assembly of quantum dots. 

As a computational expedient, each dot within an as¬ 
sembly can be treated as a lattice site, n, within discrete 
meso-scale Hamiltonian with a subscript:^^’^^ 


= ^ (7) 

n m^n 


Here Jmn is the coupling between the mth and nth sites. 
The character of exciton transport can be then delin¬ 
eated by the size of these coupling parameters relative to 
the associated reorganization energies. Suppose all of the 
dots and and their couplings to be identical so that the 
assembly is characterized by single A and J values. Then 
A ^ J leads to the the coherent transport dynamics of 
the Redfield equation^^’^^ and A ^ J results in the in¬ 
coherent, hopping transport of the Forster equation^^’^^. 
The intermediate case, wherein A ^ J, is the regime of 
partially coherent transport of interest here.^^ 

There exist two primary design challenges in creating 
quantum dots assemblies that exhibit partially coherent 
transport. The first is to reduce the reorganization en¬ 
ergy, A, as much as possible, taken on in the present 
investigation. The design of bridge assemblies with suffi¬ 
ciently high excitonic coupling, J, will be the subject of 
a future work. 
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The reorganization energy of quantum dots can be re¬ 
duced by either disrupting or constraining the exciton- 
phonon interaction because this mediates the structural 
response to electronic excitation. Each idea is explored 
with two distinct approaches intended to demonstrate 
how such mitigation can be designed into quantum dot 
solids with the aim of enabling partially coherent trans¬ 
port. 

It is not sufficient to simply reduce the reorganization 
energy, though, as the shape of the spectral density, 3, 
has a strong influence on the rate, i?, at which the ex- 
citonic superposition loses phase coherence. This can be 
elucidated by considering a SiQD dimer,for which 
this decoherence rate can be expressed as sum of a pure 
dephasing, Rd, and a relaxation, R^: 

R = Rr T Rd 

Rr - {2jfS{h)/{2h'^) (8) 

Rd = fe25(0)/(262), (9) 

Here iP' = 8e^ + (2 J)^, 5e is the energy difference between 
two SiQDs, S{uj) = 3(<^)coth(fici;/(2/cBT)) is the ther¬ 
mally weighted spectral density, and is Boltzmann’s 
constant. The dephasing expression determines the rate 
of decay of off-diagonal elements of the density operator 
after tracing over all of the phonon states. In contrast, 
the relaxation expression corresponds to the population 
transfer between eigenstates of the relevant system that 
are accompanied by energy dissipation into the phonon 
reservoir. 

If the two QDs are identical in this dimer system the 
pure dephasing rate Eq. 9 would be zero since fe = 0 
so that the decoherence rate is determined solely by the 
relaxation rate of Eq. 8. Within this setting, the rate of 
decoherence is proportional to the spectral density at a 
vibrational energy equal to the excitonic coupling, J. It 
is therefore important to suppress, as much as possible, 
the spectral density for vibrations in this range. 

Eor nonzero fe, on the other hand, the decoherence 
rate is the sum of rate Eq. 9 and Eq. 8. In that case, 
the pure dephasing rate is proportional to the slope of 
the spectral density.This implies that a spectral den¬ 
sity with a dearth of exciton-phonon coupling at the low- 
energy end of the vibrational spectrum will result in a 
flat character there and a small slope at the origin. Both 
of these key characteristics of the spectral density are 
illustrated in Eigure 3. 


II. COMPUTATIONAL APPROACH 

All ground state structures were optimized using 
density functional theory (DET) as implemented in 
DMOL.^^ An all-electron approach was used with ex¬ 
change and correlation effects accounted for by the gener¬ 
alized gradient approximation (GGA) parameterized by 
Perdew, Burke and Ernzerhof (PBE)^^ for SiQDs and 
hybrid B3LYP^^ for organic molecules. A real-space. 



FIG. 3. Schematic of optimized shape of spectral density, 3, 
in which there exists a small slope for low-energy vibrations 
and small reorganization energies for vibrational modes in 
the vicinity of the dot-to-dot electronic coupling, J. Along 
with a low reorganization energy, A, high electronic coupling, 
J, spectral densities with this character will lead to longer 
lifetimes for excitonic superpositions. 


double numeric plus polarization (DNP) basis was used 
along with an octopole expansion specify the maximum 
angular momentum function. Optimization was deemed 
converged when the energy difference between successive 
configurations was less than 2 x 10“^Ha and the maxi¬ 
mum force on any atom was less than 4 x IQ-^Ha/A. 

While TDDET and Cl are able to calculate the low¬ 
est excitation state and perform geometry relaxation to 
its equilibrium geometry, these are impractical for large 
systems since they are too expensive. Instead, a ASCE 
method was employed because it is no more computa¬ 
tionally expensive than a small set of ground state DET 
calculations. The method is based, in fact, on DET in 
which an iterative procedure is used to obtain a set of 
Kohn-Sham orbitals: 




( 10 ) 


The electron density is the trace of the density operator, 
constructed using the first N occupied kets: 


N 

P = ^\Pi){Pi\, P = Tr(p), (11) 

i=l 

In standard ASCE,^^’^^ a second iterative analysis is 
carried out but with the electron density operator, p, 
now constructed using the lowest N-1 orbitals and the 
orbital. With this change (A) imposed on the 
structure of the density, self-consistent field (SGE) it¬ 
eration is then carried out on Eq. 10 as usual, hence 
the name ASCE. Both accuracy and the rate of conver¬ 
gence can be improved, though, by performing this SGE 
iteration using a density operator comprised of the low¬ 
est N-1 orbitals and a linear combination of unoccupied 
orbitals.This is especially important when there is 
degeneracy at or near the HOMO or LUMO. Our highly 
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symmetric quantum dots typically exhibit such degen¬ 
eracies, so we use Xjyj A^deg Z^naeg I'^^deg) 

(A/'+l)^^ Kohn-Sham orbital, where A^deg is the degener¬ 
acy of corresponding state. Carrying out iteration pro¬ 
cess of Eq. (10) with this modified density results in 
the Franck-Condon excitation energy, point A in the left 
panel of Figure 2. To calculate the reorganization energy, 
a geometry optimization was subsequently performed to 
obtain point B in the left panel of Figure 2. 

This approach does not address the excited state spin, 
though, which must be antisymmetric (singlet), and some 
care must be taken to estimate this properly. We first 
assume that all N orbitals are spin restricted so that the 
promotion of a Kohn-Sham particle from orbital N to 
orbital N-\-l implies that the spin state is determined by 
particles in orbitals N-1 and These two orbitals 

can then be used to construct the requisite singlet wave 
function: 

- i;t)). (12) 

As is typical in electronic structure codes, though, it is 
not possible to assign such a superposition of spin pairs 
directly. However, the excited singlet state can also be 
written as 

iV’o) = V2m) - iV’iT), (13) 

where the triplet spin states are 

\^I) = in) 

= T(iti) + n)) 

= in) • ( 14 ) 

This implies that the desired excited state singlet energy 
surface, (green curve in both panels of Figure 2) can 
be constructed from the degenerate triplet surface, 

(red curve in right panel) and that of state |^t) (blue 
curve in right panel) 

{S\H\S)=Es = 2E^^-Et. (15) 

The energy difference between points A and B of on this 
constructed singlet surface is the reorganization energy, 
A. Since the entirety of this curve is not available in prac¬ 
tice, though, we approximate the reorganization energy 
as the difference between Ea' and Eb' on the |ti) sur¬ 
face (blue) but also calculate Ea” and on the triplet 
surface (red) for comparison. Note that an analogous ap¬ 
proach, building the desired energy surface out of those 
associated with |ti) and will not yield an energy 

relation analogous to Eq. 15 which results from one of 
the constituents being an eigenstate of the Hamiltonian. 

The spectral density, 3, of Eq. 5 is easily constructed 
once the reorganization energy. A, has been spectrally de¬ 
composed into its components, A^. To implement this, a 
vibrational (Hessian) analysis is performed on the ground 


state structures and the displacement associated with re¬ 
organization is expressed as a linear combination of the 
phonon modes. Anharmonic corrections are also required 
for the smallest dots (Si 29 H 36 ) but were found to be neg¬ 
ligibly small for all larger sizes. 


A. Estimated Accuracy of Approach 

The predictions of our computational methodology 
for calculating reorganization energy can be compared 
directly with spectral measurements by measuring the 
Stokes shift-i.e. the energy difference between absorption 
and emission peaks. The reorganization energy is then 
half of the Stokes shift. In Table I, predicted reorganiza¬ 
tion energies for common organic molecules are compared 
with such experimental measurements with nonpolar sol¬ 
vents, chosen because the solvent makes a negligible con¬ 
tribution to the reorganization energy. The average dis¬ 
crepancy between ASCF and experimental measurement 
is 84 meV. 


TABLE I. Comparison of predicted and measured reorgani¬ 
zation energies along with the magnitude of the discrepancy 
between them. Perylene is in toluene and all other molecules 
are in cyclohexane. The energy unit is eV. 



ASCF 

Exp. 

IDiff.l 

Benzene 

0.26 


0.08 

Naphthalene 

0.32 

0.30'^'‘ 

0.02 

Anthracene 

0.18 

0.2^“^ 

0.03 

Biphenyl 

0.61 

0.36®^ 

0.25 

Perylene 

0.12 

0.16®® 

0.04 


While no experimental data is yet available for the 
reorganization energy of small SiQDs, several compu¬ 
tational approaches have been used to predict val¬ 
ues as a function of dot size: Quantum Monte Carlo 
(QMC)^^, Time-Dependent Density Functional Tight- 
Binding (TDDFTB)^^ and Time-Dependent Density 
Functional Theory (TDDFT)^^. It is also common to 
use a rougher sort of ASCF approach wherein ground 
state DFT orbitals are calculated and then the excited 
state geometry/energy is obtained via a second DFT cal¬ 
culation in which the ground state HOMO is empty and 
the LUMO is filled. Table H shows the comparison of 
our calculated reorganization energy with these litera¬ 
ture results where the rougher ASCF scheme is referred 
to as DFT. With QMC taken to be the most accurate, 
the average ASCF difference from QMC is 57 meV. 
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TABLE IL Reorganization energy calculated by ASCF, DFT, 
QMC, TDDFTB and TDDFT. The energy unit is eV. 



ASCF QMC®® 

TDDFT®® 

TDDFTB®’^ 

DFT®® 

SbgHae 

0.49 

0.50 

0.82 

0.92 

0.35 

SissHae 

0.27 

0.40 

0.72 

0.74 

0.29 

Si66ll64 

0.16 



0.52 

0.25 

Si29H24 

0.23 

0.20 



0.17 


III. DESIGNING FOR REDUCED 
REORGANIZATION ENERGY 


The reorganization energy of isolated quantum dots 
can be reduced in one of two ways: by minimizing the 
affect of exciton creation on structural re-arrangement; 
and/or by confining the dot so as to inhibit the restruc¬ 
turing that follows excitation. Four design strategies 
were explored that focus on one or the other of these 
primary objectives. Mitigation of the affect of the ex¬ 
citon on bonding was explored by: (Al) increasing dot 
size; and (A2) functionalizing the dot surface so as to 
modify the frontier orbitals. The influence of structural 
confinement on reorganization energy was quantified by: 
(Bl) surface reconstruction of the hydrogen bonding; and 
(B2) encapsulation of the dot within an oxide shell. 

The most straightforward way of reducing the exciton- 
phonon coupling is by spreading the exciton out over 
more atoms. In polyacenes, for instance, the excitonic 
reorganization energy decreases monotonically with in¬ 
creasing number of phenyl rings.This is due to the ex¬ 
tended TT conjugation of such molecules, but distributed 
orbitals typify defect-free, crystalline quantum dots as 
well. Our computational results quantify an analogous 
trend with increasing dot size, summarized in Figure 4 
and Table III. It was further found that there exists an 
inverse linear dependence between excitonic reorganiza¬ 
tion energy and the number of Si atoms as shown in Fig¬ 
ure 5. This is consistent with a relationship previously 
noted for charge transfer excitons in CdSe dots^^. As 
noted in the discussion of computational methods, the 
reorganization energy as the difference between Ea' and 
Eb> on the |ti) surface (blue) of the right panel of Fig. 
2. For comparison, the analogous values associated with 
and from the triplet surface (red) are: 747 meV 
(Si 29 H 36 ), 317 meV (SissHse), 296 meV (Si66H64) and 82 
meV (Si 78 H 64 ). From these results it can be concluded 
that the spin state does not make a meaningful change 
to the reorganization energy for dots 1.4 nm in diameter 
or larger. 

Since having a small dot size is important for bet¬ 
ter transport, fewer defects, and better control of op¬ 
tical properties, consideration was next given as to how 
to reduce exciton-phonon coupling for dots of fixed size. 
This issue has been considered for photosynthetic organic 
molecules where reorganization energy can be reduced by 



FIG. 4. The reorganization energy of H-terminated SiQDs 
changing with number of Si atoms. 
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FIG. 5. The linear dependence of reorganization energy on the 
inverse of the number of Si atoms per quantum dot, Nsi- An 
equivalent relation is also given in terms of the dot diameter, 
F), measured in nm using the bulk diamond Si volume per 
atom of 0.02 nm^. 


minimizing the role that frontier orbitals play in bond¬ 
ing. This can be engineered through the addition of elec¬ 
tron withdrawing moieties that contribute to the frontier 
orbitals without significant modification of the bonding 
orbitals. As a proof of concept, a phenalenyl radical with 
such properties was found to have record low reorgani¬ 
zation energy among organic molecules of similar size.^^ 
Such character also underlies the efficiency of lanthanide 
nanoparticles in carrying out excitonic energy transfer 
up conversions.^^ There the partially filled 4f level plays 
the role of frontier orbitals while bonding is dominated 
by more spatially distributed, filled 5s and 5p states. 
This strategy can be applied to SiQDs by replacing H- 
termination with an electron-withdrawing group. Sim¬ 
ple Cl termination was deemed sufficient to illustrate the 
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idea, and Figure 6 shows that, indeed, this functional¬ 
ization mitigates what had been a strong bonding role 
played by the HOMO. The LUMO exhibits a similar 
shift. For the smallest dots, Si 29 H 36 , a drop in reorga¬ 
nization from 494 meV to 67 meV resulted. The results 
for all dot sizes, provided in Table III and plotted (blue) 
in Figure 9, indicates that there is not a strong size de¬ 
pendence to this effect. 



1.5 nm 




1.5 nm 



FIG. 6. (a) HOMO distribution of H-terminated Sirs- (b) 

HOMO distribution of Cl-terminated Sirs. 

With two methods explored for reducing the exciton- 
phonon coupling, we now turn attention to the second 
design strategy for reducing excitonic reorganization; the 
application of constraints that prevent or minimize the 
restructuring despite a driving force to do so. This too 
has precedence in organic systems. Fuller cues, for in¬ 
stance, have an excitonic reorganization energy that is 
particularly low because they are extremely rigid.This 
is one of the reasons that they figure so prominently 
in organic photovoltaic materials such as P3HT/PCBM 
blends.The excitonic reorganization energy of phthalo- 
cyanines, already low, can be further decreased with 
dimerizing end groups that increase rigidity. 

This concept is explored for SiQDs with two types of 
constraint engineering. The first exploits the fact that 
dihydride termination of surface silicon atoms can be re¬ 
duced to monohydride bonding by decreasing the partial 
pressure of environmental hydrogen.This results in ad¬ 
ditional Si-Si bonding (Figure 7) that has the structural 
character of a thin, stiff skin. Such surface reconstruction 
reduces the reorganization energy as detailed in Table HI 
and plotted (green) in Figure 9. Note that data was not 
generated for SissHse since the H atoms are not suffi¬ 
ciently close together for surface reconstruction to occur. 

A second approach for constraining structural rear¬ 
rangement is to encapsulate SiQDs within a rigid matrix. 
To test this idea, an Siss dot was encapsulated within a 
monolayer of Si02. As shown in the center of Figure 
8, this reduces the excitonic reorganization energy from 
271 meV to 161 meV. To develop a sense for the impact 
a thicker oxide shell would make, we also calculated the 
reorganization energy for SissHse with the H atoms held 
fixed. This resulted in a remarkably low excitonic reor¬ 
ganization energy of 36 meV, as shown at right in Figure 
8 . 


FIG. 7. (a) Si 78 H 64 and (b) the same quantum dot after sur¬ 
face restructuring produces 6 H 2 molecules leaving Si 78 H 52 . 



FIG. 8. The reorganization trend as Sias becomes more 
rigid. The blue dots corresponds to H-terminated Siss; the 
green dot corresponds to Sias within one oxidized Si shell; 
the red dot corresponds to H-terminated Siss with H fixed to 
mimic Siss within an infinite large matrix. 


A summary of all four analyses is provided in Table HI 
and Figure 9. The lowest reorganization energy achieved 
is 17 meV (137 cm“^) for a surface reconstructed SiQD 
with a diameter of 1.7 nm. 


TABLE III. Reorganization energy of H-terminated, Gl- 
terminated and surface-restructured SiQDs. The energy unit 
is meV. 


SiQD diameter(nm) 

0.9 

1.1 

1.2 

1.4 

1.7 

# Si Atoms 

29 

35 

66 

78 

147 

H-terminated 

494 

271 

164 

80 

40 

Gl-terminated 

67 

57 

98 

46 

25 

Restructured 

226 

- 

83 

18 

17 


We now turn to a consideration of the shape of the 
spectral density to determine the extent to which it can 
be optimized along the lines of Figure 3. As a specific 
example, the spectral density of Si 78 H 64 dots is compared 
with that of its restructured counterpart, Si 78 H 52 . Figure 
10 shows both spectra along with their Drude-Lorentz fits 
using Eq. 6 with 34 peaks for Si 78 H 64 and 13 peaks for 
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FIG. 9. Plots of excitonic reorganization energy versus SiQD 
diameter for: standard H-termination (red); surface recon¬ 
structed dots with lower H content (green); and electron- 
withdrawing Cl-termination (blue). The connecting lines are 
a guide to the eye. 


Si 78 H 52 . Reconstruction modifies the lower energy end 
of the spectrum in a way that should significantly reduce 
the decoherence rate, and restructuring also suppresses 
the spectrum in the physically desirable 500-700 cm“^ 
range of electronic couplings. This too would lead to 
a system more capable of supporting partially coherent 
transport. 


IV. CONCLUSIONS 

Sufficiently low dielectric screening in quantum dot as¬ 
semblies makes charge dissociation improbable following 
photo excitation. Energy is then transported via exci- 
tons in what is typically a random walk. Such diffusive 
propagation is the result of a loss of coherence in the 
superposition of excitonic states that are present in the 
initial photo-absorption, but maintaining an element of 
such coherence is thought to result in a robust, more effi¬ 
cient form of energy transport. This requires a low reor¬ 
ganization energy for isolated dots and a large excitonic 
coupling between dots. The current work has focused on 
ways in which reorganization energy can be reduced by 
either reducing or constraining the associated exciton- 
phonon coupling. Two design paradigms were explored 
for each of these. 

Reduction in the coupling itself can be accomplished 
most simply by spreading out the electronic excitation 
over more atoms, for instance by increasing dot size. The 
reorganization energy. A, was found to be inversely pro¬ 
portional to the cube of the dot diameter, D, and a simple 
relationship was identified: A = 236/D^ with D in nm 
and A in meV. Since having a small dot size is important 
for better transport, fewer defects and easier functional¬ 
ization, though, a second way of reducing exciton-phonon 
coupling was also explored. Dots were functionalized so 



FIG. 10. Spectral density of Si 78 H 64 (red) compared with 
that of its surface-restructured counterpart, Si 78 H 52 (green). 
Vertical lines indicate non-negligible contributions to the re¬ 
organization energy and have been normalized to the peak 
value of the fitted spectral density. 


that the frontier orbitals, of which excitons are primar¬ 
ily comprised, are not the bonding orbitals that generate 
phonons. This idea has been successfully demonstrated 
for small organic molecules, and it is also underlies the ef¬ 
ficiency of lanthanide nanoparticles in efficiently carrying 
out excitonic energy transfer upconversions.^^ Our anal¬ 
ysis indicates that this is also possible for SiQDs using Cl 
atoms-not because Cl is considered a candidate function¬ 
alization but to demonstrate proof of the concept. For 
the smallest dots, a drop in reorganization from 494 meV 
to 67 meV resulted. This is a surface effect, however, and 
further reduction in reorganization energy is minimal as 
dot size is increased. 

The second design paradigm is to constrain the system 
so that it cannot carry out the structural re-arrangements 
that generate reorganization energy. This was first ex¬ 
plored by controlling the hydrogenation of dot surfaces 
which can be carried out by changing the partial pressure 
of hydrogen in the environment.^^ The results indicate 
the reorganization energy is significantly reduced as a re¬ 
sult of what amounts to a stiff skin on the surface of the 






















dot that resists the driving force to reorganize. As with 
the standard H-terminated dots, this reorganization en¬ 
ergy decreases with increasing dot size, and 17 meV was 
the minimum value computationally measured. Matrix 
encapsulation is also a promising avenue for reducing the 
motion that generates the reorganization energy well. A 
monolayer of Si02 around Siss reduced A from 271 meV 
to 161 meV and freezing the H-atoms on the dot surface 
dropped this value down to only 36 meV. 

This computational study of 1-2 nm silicon quantum 
dots is not focused on reducing excitonic reorganization 
energy as much as possible. Rather, it identifies ways in 
which this energy can be reduced in order to offer direc¬ 
tion to future work. Even so, the results suggest that it 
is reasonable to identify a goal of creating 1-2 nm silicon 
quantum dots with reorganization values in the range of 
1-10 meV. This overlays the range typically considered 
in association with photosynthetic complexes of 0.1 to 
10 meV^^-naturally occurring assemblies that serve as 


a bio-inspiration for what may be improved upon with 
engineered organic/inorganic nanostructures.^^“^^ More 
important than the range itself, though, successful de¬ 
signs will need to produce excitonic coupling of the same 
order or larger to create quantum dot solids with a mea¬ 
sure of wave-like character. 
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